The charge and low-frequency response of normal-superconducting heterostructures 
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The charge distribution is a basic aspect of electrical transport. In this work we investigate the self- 
consistent charge response of normal-superconducting heterostructures. Of interest is the variation 
of the charge density due to voltage changes at contacts and due to changes in the potential. We 
present response functions in terms of functional derivatives of the scattering matrix. We discuss 
corrections to the Lindhard function due to the proximity of the superconductor. We use these results 
to find the dynamic conductance matrix to lowest order in frequency. We illustrate similarities and 
differences between normal systems and heterostructures for specific examples like a ballistic wire, 
a resonator, and a quantum point contact. 
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I. INTRODUCTION 

During the past decade mesoscopic systems consisting 
of both normal and superconducting parts have attracted 
considerable attention. Microscopically, the interesting 
physics stems from Andreev reflection. An incident par- 
ticle is reflected as a hole and a Cooper pair is gener- 
ated in the superconductor. This results in an effective 
charge transfer of 2e and correlations between Andreev 
reflected electron hole pairs (the proximity effect). These 
effects have been investigated in many experimental and 
theoretical works |l|,g| focusing mainly on the station- 
ary transport regime (dc-conductance) [0,Ef and the low- 
frequency noise (shot noise) la, ok The ac-regime has at- 
tracted much less attention fK- HI . 

In an Andreev process the electron and hole parts of 
the wave function contribute with opposite charge. It is 
therefore interesting to investigate the low-frequency ac- 
transport of NS-systems, since this problem requires a 
electrically self-consistent discussion of the charge distri- 
bution in the sample. This self-consistency is of impor- 
tance not only for ac-transport but also for the discus- 
sion of charge fluctuations and the non-linear transport 
regime pT[ . 

In this work we have in mind the interplay of two main 
properties of hybrid structures: On one hand raising or 
lowering the voltage at a normal contact of the sample 
will not inject an additional charge into regions where the 
wave functions contain electron and hole amplitudes of 
equal magnitude. This is in strong contrast to a purely 
normal conductor! On the other hand screening is a prop- 
erty not only of the states at the Fermi surface but of the 
entire electron gas. Thus the ability of a hybrid structure 
to screen an additional charge is essentially the same as 
that of a normal conductor. 

Our results show two main differences compared to 
purely normal systems: first, the coupling of carriers 
with opposite charge reduces the interaction with nearby 
gates. Second, Andreev reflection increases the dwell 
time inside the structure and this affects the ac-response. 

The paper is organized as follows: In section Q|, we 



derive an expression of the charge density in terms of 
functional derivatives of the scattering matrix. We next 
discuss the charge density response to external and in- 
ternal potential perturbations. In section IV we use 
these results to formulate a self-consistent theory of low- 
frequency ac-response. To illustrate our results we con- 
sider in section M a number of examples: a ballistic wire, 
a resonant structure and the quantum point contact con- 
nected to a superconductor. 



II. CHARGE DENSITY 




FIG. 1. A mesoscopic scattering region is attached to N 
normal reservoirs and one superconducting reservoir. It is 
described by two scattering matrices sn describing the normal 
conducting region and sa describing the Andreev reflection. 
Each terminal has its own bias voltage V a . The electrostatic 
potential U inside the scatterer is calculated self-consistently. 



In this section we derive the general expression for the 
charge density (p(r)) in the scattering problem sketched 



in Fig. W. A scattering region is attached to N normal 
leads and one superconducting lead. Every normal lead 
is characterized by its applied voltage V a , the supercon- 
ducting lead by its pair potential A and the bias Vs- For 
all the calculations we may choose Vs = 0. The fact that 
we only allow for one superconducting terminal excludes 
all time-dependent Josephson-like effects. For an intro- 
duction to the applied formalism we refer the reader to 
Ref. p3] . The whole system is described by its scattering 
matrix 
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The element sJL for example is the current amplitude of 
a hole that leaves through lead a and has entered with 
unit current amplitude as a particle through lead /3. We 
represent each scattering channel by its own lead to save 
two indices. 

It is conceptually useful |13| to imagine the scattering 
matrix s being assembled from a part that describes the 
reflection and transmission in the normal region s p Jj and 
from a part that describes the Andreev processes at the 
interface between normal metal and superconductor sa- 
Only this second matrix leads to coupling between the 
particle and hole scattering states. However, the follow- 
ing derivations do not depend on this assumption. For 
energies E smaller than the gap energy | A| only reflection 
takes place at the interface. The total scattering matrix 
s then has the dimension (2Nn) 2 where A^v = J2 N a 
is the total number of channels leading to normal reser- 
voirs. Above the gap we must also include transmis- 
sion processes, and therefore the dimension changes to 
(2N + 2 As) 2 . The superconducting terminal adds iVg 
more channels. 

For the following it is helpful to introduce local partial 
densities of states (LPDOS) gl 
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This expression is valid for one channel per lead. A true 
multichannel expression would include a trace over the 
channels. The value u(l h ,r p ,2 p ) for example describes 
the density of particles at location r that entered as par- 
ticles through contact 2 and leave as holes through lead 
1. In definition (0) we denote the quasi-particle charge 
by q p ' h = ±q. The LPDOS must be calculated as func- 
tional derivatives of the scattering matrix with respect 
to the electrostatic potential U. To gain the information 
about particles and holes separately |l4J we artificially 
split up the electrostatic potential U in a part that acts 
on particles U p and another that addresses holes U h . The 
Bogoliubov-de Gennes Hamiltonian then takes the form 



H = 



H +qU p (x)-E F A 

A* -(H Q +qU h (x)-E F y 



(3) 



These equations have to be solved including a small vari- 
ation of the electrostatic potentials U p and U h in order to 
get the scattering matrix and its functional derivatives. 
The above defined LPDOS are not independent. On 
one hand, the LPDOS obey reciprocity relations. This 
has been investigated in reference |1J] . On the other hand 
the particle-hole symmetry of the Bogoliubov-de Gennes 
equation implies 



u(a v ,r K ,/3^,E) = v(ap,r R ,Pp.,-E). 



(4) 



The bar denotes the opposite (p/h — h/p). Both symme- 
tries can be used to reduce the expense of the calculation. 
The charge density inside the normal-superconducting 
heterostructure can be entirely expressed by the LPDOS 
(and therefore by the scattering matrix) and the occupa- 
tion factors of the attached reservoirs 



(p(r)) = ljr oo dEj: a 
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(5) 



The occupation factors include the bias voltage of the 
normal terminals /£' (E) — f(E — q p ! h Vp). Here / is 
the Fermi function. Note that the occupation factors 
vary in opposite directions for particles and holes. In Eq. 
(H) we have double counted the particle-hole excitations 
and hence drop a factor of two for spin degeneracy. The 
derivation of this result is outlined in appendix [A|. 



III. CHARGE RESPONSE AND GAUGE 
INVARIANCE 

Given formula (m we are now in a position to calculate 
the charge density response 5p(r) to both internal and 
external potential variations 



Sp(r) 
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dp(r) 



8Vp- / dr'fl(r,r')SU{r'). (6) 



The first contribution is the bare charge injected from 
the leads due to the shift of the occupation factors, and 
is proportional to the injectivities dp(r)/dVp. The sec- 
ond contribution arises from the change of the internal 
potential due to screening (the potential itself will be de- 
termined in the following subsection), and involves the 
Lindhard function ft (r, r') = —dp(r)/dU(r'). 

The injectivity from the normal leads can be calculated 
straightforwardly from the charge density dq) 






g VK,r K ,^) (7) 



and depends at low temperatures as expected only on 
properties at the Fermi energy. The other quantities con- 
tained in the balance equation (0) need a more careful 



analysis. Its technical details are explained in appendix 



IV. LINEAR RESPONSE CALCULATION 



The procedure of calculating the non-local Lindhard 
function n(r, r') leads to second order functional deriva- 
tives that cannot be simplified further. However, if we 
assume the Lindhard function to be local, II(r, r') — 
S(r — r')lUr), we can express it by the above calculated 
LPDOS (g). This assumption is correct if the electro- 
static potential varies only slowly on the scale of the 
Fermi wavelength Xp. To express the result in a transpar- 
ent way we split up the energy E in the same way as the 
potential U(r). We introduce the quantities E p l h that 
denote the energy dependence of the scattering matrices 
for particles and holes s^v (E p / h ). Furthermore we need 
a critical energy E c with the property |A| <C E c <C Ep. 
Such an energy is useful, since particles and holes with 
energies outside the range [—E c , E c ] have a negligible An- 
dreev reflection probability that decays as ~ \A\ 2 /E 2 . In 
a short calculation given in appendix we find for the 
local Lindhard function 

n ( r ) = q 2 Y, a al l '( a P' r p>Pp>- E c) + v( c Xh,rh,Ph,+Ec)] 
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If the local electrostatic potential U(r) fulfills the condi- 
tion Ep — U(r) 3> | A | the first line of Eq. (||) dominates. 
Therein we can replace the sum over the LPDOS by the 
local density of states at E — for the equivalent nor- 
mal conducting structure. This is the expression for the 
Thomas-Fermi screening in a purely normal sample. The 
screening properties are not affected by the presence of 
the superconductor. Our result becomes more interesting 
for Ep — U(r) ~ |A| when the second and third term of 
Eq. rt|) contribute significantly. We later illustrate this 
case with a specific example. 

A simple argument allows us to get the injectivity from 
the superconducting terminal dp(r)/dVs without any 
further calculation. The Bogoliubov-de Gennes equa- 
tions are gauge invariant, a simultaneous change of all 
external and internal potentials by the same amount will 
not lead to any charge inside the system. Setting the left 
side of Eq. (ra) to zero gives therefore 



H(r) 



0=S,1..N 



dp(r) 
dVfi ' 



Since IL(r) and the injectivities of all normal contacts 
are known we can use this relation to find the injectivity 
dp(r)/dVs of the superconducting contact. 



In order to get the low-frequency ac-response of our 
system it is necessary to distinguish two contributions to 
the current. On one hand we have the current flow I bare 
of non interacting particles which can be accessed by a 
linear response theory. On the other hand we may not 
neglect the screening currents I scr due to interactions. 
The low-frequency conductance matrix can be generally 
written as 
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(10) 



where the "emittance" matrix E consists of two parts 

77 rpbare _i_ rpscr 

The screening currents may be calculated quasi-stati- 
cally solving a Poisson equation self-consiste ntly . This 
procedure is described in detail in reference |15f| . Here 
we cite only the result 



rpscr 
E a /3 



drdr' g(r, r') 



dp(r) dp(r') 
~dK dVa ' 



(11) 



which is valid in the presence of time reversal symmetry. 
The kernel g(r, r') is given by J dr"AnIi(r, r")g(r",r') — 
W 2 g(r,r') = An8(r — r'). In a discretized model the 
Laplace-operator may be replaced by a capacitance ma- 
trix. 

To find the bare contribution E bare to the emittance we 
proceed as in reference [|16| . We use the current operator 
at the normal conducting terminal a 



J = IV « f dEdE'e +i{ - E ~ E '} t / h 



(12) 



a^(E)a%(E')xAffi(a,v,E,E'), 



in a simplified form valid in the low- frequency limit. 
The full current operator has been given for example 
in reference |l7]. In Eq. ( |l2] ) the indices a,P,/3' de- 
note leads (and channels), v,w,w' distinguish particle 
and hole states. The operator d\{E) for example cre- 
ates a hole of energy E incident into lead 1 . The matrix 
elements A are given by 



Appi (a v ,E,E ) = 5 a p5 a p'5 V w8v 
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(9) The ac-conductivity can then be obtained from 



Qaa' 



T)lU 
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/a(T),Ia'(0) 
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The evaluation of the commutator is mostly straight for- 
ward. As in previous works |16| we use the unitarity of 
the scattering matrix and the thermal occupation of the 
reservoirs. As in the case of a purely normal system we 
are left with a doubled energy integral. We can evaluate 



this integral through a path deformation in the upper 
complex plane where the scattering matrix is analytical. 
In the end we expand the result up to first order in fre- 
quency. The result for the dc-conductance 
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is identical to the one established in the literature 
0,|8Jl3[l9]. This serves as a check of our calculation. 
In Eq. (|15[), T^ t v „ = (s^i)* s^~ is the transmission proba- 
bility from channel (3 to channel a. The bare emittance 
can be expressed by global partial densities of states (GP- 
DOS) 



aS a/3 
8E ' 



(16) 



and becomes 

E b a y = (dE (-|£) EaW^'/U (i7) 

This equation shows that the bare emittance may change 
its sign. This simple calculation provides only the emit- 
tance matrix elements between normal terminals. A di- 
rect calculation of the current at the superconducting 
reservoir would involve a self-consistent evaluation of the 
pair potential in the superconductor. Its phase "carries" 
the supercurrent. Nevertheless, the missing elements of 
the emittance matrix can be reconstructed from the con- 
ditions 



2_, Ea P -/2 Ea 











(18) 







that express current and charge conservation. 



V. EXAMPLES 

We now present some simple calculations to illustrate 
how the presence of a superconducting terminal affects 
the ac-properties of a mesoscopic sample. We emphasize 
that these examples are not designed to model a realistic 
sample completely, but should exhibit qualitatively the 
main features. 



A. Ballistic Wire 

As a first example we discuss briefly the emittance of a 
ballistic wire with one open channel at zero temperature. 
The results can be easily generalized to more than one 
channel. The geometry of the sample is shown in Fig. g. 
The wire is attached to two reservoirs (1,2). Reservoir 1 
is always normal conducting. The second may be normal 
conducting, superconducting or completely disconnected 
from the wire for comparison. 



a) 



b) 




FIG. 2. Examples: a) A one channel ballistic wire is at- 
tached to a superconductor. Its length is denoted by £, the 
internal interaction modeled by a capacitance per unit length 
c between wire and gate. The emittance is enhanced by a 
factor of 4 compared to a purely normal wire, b) The same 
geometry for a quantum point contact which is described by 
an electrostatic dipole of capacitance Co. In the low trans- 
parency limit the ac-conductance is dominated by the geo- 
metric capacitance Co- In the opposite limit the quantum 
point contact shows enhanced inductive behaviour. 



The wire is described by its length £ and the dimen- 
sionless parameter a = 4q 2 /hvp which describes the DOS 
at the Fermi level in the wire. Note that a is also the lo- 
cal Lindhard function of Eq. (ph. The interaction in the 
wire is modeled by a third gate terminal (3). It is coupled 
to the wire by a geometrical capacitance per unit length 
c and assumed to be macroscopic. Thus we can replace 
the Laplace operator in Eq. (O) by —inc. For a detailed 
description of this system see for example pfj . 



En 

-El 2 

E22 



normal conducting superconducting disconnected 



4(c+a) 

al(a + 2c) 
4(c+o) 

a 2 t 
4(c+a) 



-a(t + to/2) 

a(£ + to/2) 



W-«5o/2 



(c+a) 






TABLE I. Comparison of the emittance elements of a bal- 
listic wire connected to two normal reservoirs, a normal and 
a superconducting reservoir and to one normal reservoir only. 



As a last parameter we need the coherence length of 
the superconductor £ = ftuj?/|A|. We neglect the self- 
consistency of its pair potential. 

Tabic | summarizes the results for the three cases. The 
missing elements of the emittance matrix can be recon- 
structed from Eq. (jig). We add some observations to 
explain the differences between the results. The response 
of the disconnected wire is purely capacitive, while the 
open wires act inductively. In the limit of charge neu- 
trality c <C a the inductive emittance of an open wire 
En grows by a factor of four in the presence of a super- 
conductor. On one hand the bare emittance is doubled, 
because an incoming electron leaving as an Andreev re- 
flected hole stays twice as long in the wire. On the other 
hand this effect is not weakened by a capacitive screened 
emittance, because the injectivity from the normal lead 
into the wire dp/dVi is zero. This leads to another factor 
of two. Additionally, the evanescent quasi-particle wave 
contributes to the bare emittance, the wire acquires an 
effective length £ + £ /2 (we use the assumption that the 
Fermi velocities are the same on both sides of the NS- 
interface) . 

The emittance E13 = —En — E12 is always zero in the 
presence of a superconductor. The gate and the normal 
terminal are only connected via the capacitance. But 
this capacitance cannot be charged from the normal side 
1 because the above mentioned injectivity dp/dVi is zero. 
Every injected electron comes back as a hole that com- 
pensates its charge. Therefore, the ac-response to a bias 
at the normal end must be zero. Vice versa the capaci- 
tive element £23 becomes twice as big because of a dou- 
bled injectivity dp/ 'dV<i in the limit of potential neutrality 
O a. 



B. Resonator 

As briefly discussed after Eq. (||) we expect the in- 
ternal charge response to be very different from a purely 
normal system, if the states inside the superconducting 
gap give the dominant contribution to the Lindhard func- 
tion. To further investigate this point we now consider 
the charge on a series of resonances of Breit-Wigner type 
coupled to a superconductor. For analytical simplicity 
we take the resonances equidistant. A single resonance 
of this kind is discussed for example in pll . Its exper- 
imental equivalent might be a level on a quantum dot 
which is coupled to one normal and one superconducting 
lead through tunneling barriers. However, this situation 
is normally treated in a Coulomb blockade picture where 
the strongly fluctuating electrostatic potential is treated 
as an operator. In our theory the potential becomes a 
c-number. Therefore, this model is limited to the case of 
a small charging energy. 

The resonances are described by the following scatter- 
ing matrix in the normal part 



s N,ij — 01 



yr~r~ E - qU + i(T N + r s )/2 

ITT-^- COt7T — . 



(19) 



The coupling width of the resonances to the normal lead 
is Tn, to the superconducting lead it is Ts- The level 
spacing on the resonator is given by Al . The above given 
scattering matrix is only unitary in the limit T N ig <C A^ 
which means for small coupling. 




El. Potential qU / A L 



FIG. 3. A resonator with level spacing Al is attached to 
a normal and a superconducting wire with resonance widths 
Fjv/s respectively. Their average is called T — (Ts + Fiv)/2, 
the parameter A = (r^ — Fjv)/(rs + Tm) describes their re- 
lation. The big plot shows the charge Q on the resonator 
as a function of its electrostatic potential U for different pa- 
rameters A. The lineshape changes (see Eq. (Ell)), and for 
increasing A the charging decreases. This is illustrated in the 
inset that shows the charge per resonance Qres as a function 
of A. 



This scattering matrix is periodic in the energy E and 
the integrals in the equation for the charge density (0) 
diverge. We therefore introduce a cutoff E m i n by hand 
that defines the lowest level on the resonator. This cutoff 
lies still inside the energy range of the superconducting 
gap I -Emm I <§; |A|. This allows us to use the scattering 
matrix 



sa 



„ph 



-i 

-i 



(20) 



to describe the Andreev reflection. This example is now 
constructed in such a way that the correction of the Lind- 
hard function due to the presence of the superconductor 
(second term in Eq. (pf)) becomes very important. Fig. 
g shows a numerical evaluation of the charge on the res- 
onator Q as a function of the electrostatic potential U on 
the resonator. We describe the coupling by a parameter 
A = (r s - rAr)/(r s + T N ). This parameter gets 1 if the 
resonance is only connected to the superconductor and 



— 1 in the opposite case. The inset if Fig. |j| illustrates the 
charge stored on one resonance Qres a s a function of A. 
This charge is 2q if the resonator is charged only from the 
normal side (we assume spin degeneracy, this is the point 
where the Coulomb blockade gives a completely different 
picture). If the resonator is only connected to the su- 
perconductor no charge can be transferred at all. Below 
the superconducting gap |A| the states in the resonator 
are equally weighted superpositions of particle and hole 
states that cannot contribute to the charging of the res- 
onator. 

It is worth to note, that the shape of the steps in Fig. 
changes as a function of A. In the two different limits 
we obtain 



Q N = -grfarctanf^- T N ^>T S , 
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2qU 
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(21) 



N- 



These results are valid at zero temperature. For kT > 
Tn/s the steps are smeared out by the Fermi function. 



C. Quantum Point Contact 

The low-frequency conductance of a quantum point 
contact (QPC) connecting two normal leads has been 
studied in Ref. [Q. We adapt this procedure to our 
situation sketched in Fig. S. In a first step we only con- 
sider one transmission channel. We assume the QPC is 
described by a symmetric equilibrium potential. At equi- 
librium the only asymmetry stems from the presence of 
the superconducting lead. Polarization of the QPC due 
to an applied voltage leads to a dipole. (Charging vis-avis 
the gates is neglected. See however Ref. |P3| ). The size 
of this dipole is described by one single capacitance Cq. 
Furthermore, we limit ourselves to a semiclassical treat- 
ment which essentially means that the confining potential 
is sufficiently expanded in space. In this limit the second 
and third part of Eq. (||) are negligible. As a second pa- 
rameter we need the total density of states at the Fermi 
level (over a region in which the charge is not screened 
fully), when the system is entirely normal D N — q 2 Np. 
In addition, scattering at the QPC is characterized by its 
transmission probability T and its reflection probability 
R = 1 — T. As shown in [£2| the electrochemical capac- 
itance and the emittance in the purely normal system 
are 
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This result uses the fact that the semiclassical injectivi- 
ties may be written as 
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(1-R). (23) 



For example, the response dQi JdV\ originates from all 
the right going electrons in region 1 plus the left going 
ones that have been reflected at the barrier. 

It is clear that this picture will change drastically 
in the presence of a superconductor. We denote by 
Rn = 4i?/(l + R) 2 the probability that an electron is 
scattered back as an electron. The probability for An- 
dreev reflection we call Ra = 1 - Rn = T 2 /(l + R) 2 . 
The dc-conductance is of course ARAq 2 /h. The injectiv- 
ities now turn out to be 



ggf _ D 
dVi 



N 



-Rn, 



dQl 
dVi 



= 0, 



(24) 
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dV 2 
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For example in dQf/dVi, we recognize that only the elec- 
trons that return as electrons contribute to the injectiv- 
ity. We see also that the normal terminal cannot inject 
charge into the right side of the QPC which is also clear 
from intuition. 

Now we use these ingredients to find the capacitance 
and emittance of the whole QPC. For simplicity, we cite 
the results without the length renormalization due to a 
finite £o in the superconductor. We find 



C b =R 
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(27) 



In the low transparency limit (i? ~ 1) the result is the 
same as for the purely normal conducting system (J22|) . 
In the high transparency limit CR — 0) we recover the in- 
ductive behaviour of example |V A| . Again the emittance 
is increased by a factor of four in comparison to the result 

(ID- 

Fig. ^ shows a qualitative comparison of the con- 
ductance, capacitance and emittance of a multichannel 
QPC in the two geometries. We use a capacitance of 
Co = IfF and a potential U(x) = max{Vb(A 2 — a: 2 ),0} 
where A = 500nm. The constriction in y-direction al- 
lows up to five open channels with equidistant spacing 
through the contact. 

What are the restrictions of the results obtained for 
our simple model system? The assumption that the NS- 
interface is a perfect Andreev mirror seems to play the 
most important role. In this case we may neglect the 
capacitance of the NS-interface. If such a capacitance 
would be present it would decrease the inductive be- 
haviour at high transparency. 
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FIG. 4. Comparison between transport properties of a mul- 
tichannel quantum point contact either attached to a normal 
lead (narrow line) or a superconducting lead (broad line). In 
the presence of the superconductor, its dc-conductance is dou- 
bled and the low-frequency ac-response is enhanced by a fac- 
tor of four in the transparent limit. The curves are calculated 
in a WKB-approximation. The used parameters are specified 
in the text. 



shows the same capacitive behavior as its normal con- 
ducting analog in the limit of small transmission. For 
high transmission the emittance is enhanced as in the 
case of a ballistic wire. 

For the ac-conductance problem screening is necessary 
if we want to find a response that depends only on voltage 
differences and which conserves current. We have focused 
on geometries with a single NS-interface but similar con- 
siderations should apply if we deal with SNS-structures or 
more complicated geometries. Electrical self-consistency 
is relevant not only for dynamic problems but also if we 
are interested in non-linear transport or even just in the 
gate voltage dependence of stationary transport quanti- 
ties. Therefore, the considerations presented should be 
useful for a wide range of geometries and for the investi- 
gation of many different physical problems. 
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VI. CONCLUSIONS 

In this work we have extended the ac-response theory 
of normal mesoscopic conductors to hybrid normal and 
superconducting structures. This requires an investiga- 
tion of screening and a discussion of the charge density 
response to external lead voltages in the presence of An- 
dreev scattering. Global gauge invariance is valid also 
for the hybrid structures investigated here. This leads 
necessarily to the existence of an injectivity of the su- 
perconductor into the normal part of the structure. The 
charge-injectivity of the superconductor compensates the 
suppression of the charge-injectivity from a normal con- 
tact. 

Screening in hybrid structures is up to small correc- 
tions the same as in normal conductors. Nevertheless, 
the ac-response of hybrid structures exhibits marked dif- 
ferences from that of a purely normal system. For a bal- 
listic wire at one end connected to a superconducting 
reservoir, the emittance is four times as large as that of 
a purely normal wire. Furthermore, the displacement 
current induced into a nearby gate in response to an 
oscillating voltage at the normal contact (described by 
an off-diagonal capacitance element) is highly suppressed 
compared to the purely normal structure. For a resonant 
structure, the charge response to an internal variation of 
the electrostatic potential can be reduced due to the su- 
perposition of particle and hole states inside the cavity. 
A quantum point contact attached to a superconductor 



APPENDIX A: CHARGE DENSITY 

This is a short sketch of the derivation of formu- 
la (||). We express the expectation value of the 
charge density operator with help of the normalized so- 
lutions (w%(r,E))i = ((ug(r, £))*, (t#(r, £))*) of the 
Bogoliubov-de Gennes equation |0| 
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The solutions of the Bogoliubov-de Gennes equation are 
scattering states in our case. We include the usual pref- 
actors containing the group velocities in the normaliza- 
tion factors of the wavefunctions. Their mean occupa- 
tion number can be expressed by the Fermi functions of 
the reservoirs /?' (E) = f(E — <fl h Vp) depending on 
whether they describe an incoming particle p or hole h. 
qP/h _ ^ denotes their quasi-particle charge. 

The starting point of the calculation is the Bogoliubov- 
de Gennes equation (^). For the moment we allow the 
electrostatic potentials JJ v l h to be complex. The con- 
tinuity equation for the quasi-particle current j^a then 
reads 
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The complex potentials generate source terms on the 
right side of Eq. (A2). As a next step we integrate 
this equation over the volume of the scatterer. To 



this end we choose the potentials to vary like U p / h 
j_ nTP/hxf^ _ XQ y We then g e t for the current 
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In this equation we call I„ 



./out 



the total current flow into 



and out of the scattering region. Their difference is not 
zero because of the source term in Eq. ( |A2| ) . The ratio of 
both quantities can be expressed by the scattering matrix 
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The scattering matrix is a functional of the small complex 
variation ^ p /h and thus can be expanded up to first or- 
der. To evaluate the incoming current Ii n we use the nor- 



malization of the wave functions 



and get J£ in 



1/h. 



Finally we manage to express the square of the wave func- 
tions by the LPDOS given in definition (0) 



^2 / v{a v ,r p ,j3 lx ) = 
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These quantities can be inserted into Eq. ( |Al[ ) to get the 
final result (pj) given at the beginning of the article. 



APPENDIX B: LINDHARD FUNCTION 

In this appendix we explain the derivation of the local 
Lindhard function (LLF) given in Eq. (||). The non- 
local Lindhard function (NLF) is given by ±I(r, r') = 
—Sp(r)/SU(r'). We define functional potential deriva- 
tives of the LPDOS as follows 



Xaf3VKi r \) 
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Using this definition and Eq. 
as 



we can write the NLF 
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The NLF is not a Fermi surface quantity but depends 
on all energies within the conduction band. The LLF 
is a good approximation if the electrostatic potential U 
varies slowly on the scale of the Fermi wavelength. Under 
these circumstances the spatial integration appearing for 
example in Eq. (0) can be simplified 



dr'U(r,r')SU(r') = 5U(r) \ dr'U(r,r') = SU(r)U(r). 

(B3) 



To get the LLF we must therefore integrate the NLF 
over its second spatial variable r'. Because of particle- 
hole symmetry it is sufficient to keep the first part of Eq. 
( B2 ) . We may thus write the LLF in the following way 



-E e +00\ 
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where we used an energy scale E c with the property 
|A| <C E c <C Ep. Below — E c we may neglect any An- 
dreev reflection. In this case there is no coupling between 
particle and hole states and all 'crossed' LPDOS and the 
functional derivative S/5U h (r') vanish. We may further- 
more use 
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dEi 
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which holds in WKB-approximation and d/dE p = d/dE 
to get 

dr 'Xa V ff( r P' r x) - S P uS pfl S pX (--^J v {oL Vl r p ^ v ). 

(B6) 

The second integral may similarly be simplified using Eq. 
(B5). Applying partial integration with respect to the 
energy we get the result presented in Eq. <H). 
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